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At this point in time, there is very little empirical evidence on the likelihood of a space-faring 
species originating in the biosphere of a habitable world. However, there is a tension between the 
expectation that such a probability is relatively high (given our own origins on Earth), and the lack 
of any basis for believing the Solar System has ever been visited by an extraterrestrial colonization 
effort. This paper seeks to place upper limits on the probability of an interstellar civilization arising 
on a habitable planet in its stellar system, using a percolation model to simulate the progress of such 
a hypothetical civilization’s colonization efforts in the local Solar neighborhood. To be as realistic as 
possible, the actual physical positions and characteristics of all stars within 40 parsecs of the Solar 
System are used as possible colony sites in the percolation process. If an interstellar civilization is 
very likely to have such colonization programs, and they can travel over large distances, then the 
upper bound on the likelihood of such a species arising per habitable world is on the order of 10 -3 ; 
on the other hand, if civilizations are not prone to colonize their neighbors, or do not travel very 
far, then the upper limiting probability is much larger, even of order one. 


I. INTRODUCTION 

After decades of human exploration of space, it appears 
that there have been no visitors to the Solar System from 
other star systems. This is in sharp contrast to the (per¬ 
haps naive) view that, with hundreds of billions of stars in 
the Milky Way galaxy, the likelihood of finding habitable 
planets around many of these stars, and the possibility 
of life arising on a good proportion of these planets - 
as evidenced by the seeming rapidity of life developing 
on the young Earth - there should be a large number 
of technological civilizations, and at least one of these 
should have developed spacecraft capable of reaching the 
Solar System. The variance between this supposition, 
and the lack of evidence for any visitations, is known as 
the Fermi paradox HIS]. Its resolution has vexed many 
minds over the decades, although at this point, empirical 
data sufficient to address the question is lacking. 

This has not stopped the development of theoretical 
models whose intent is to provide insight into the Fermi 
paradox. An as example appropriate to our discussion, 
Landis [8] originally developed a model of interstellar col¬ 
onization processes as a simple percolation model, where 
each settled system has a given probability p of colonizing 
all of its neighboring systems; since the model was sim¬ 
ulated on a cubic lattice, this fixed the number of such 
neighbors to be six for every single star system. A gener¬ 
alization of this model was given by Hair and Hedman [5] , 
where the cubic lattice was still used, but colonization 
along diagonal connections was allowed - so each system 
had up to 26 neighboring systems - and p was the prob¬ 
ability for each neighboring system (rather than all) to 
be colonized by a neighbor. However, neither model used 
the actual “geography” of star systems; this aspect was 
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included in the model developed by Cartin @]. There a 
percolation model was also used, but instead of a fixed 
lattice, the observed positions of all known star systems 
within 40 parsecs (pc) of the Solar System were taken as 
colonization sites. 

The main focus of Ref. 4 was the following question: 
supposing each star system is the point of origin for a 
space-faring civilization with a given sociological drive 
and technological prowess, how many could reach the So¬ 
lar System with their colonization process? To codify 
this, the two parameters of the model were a probabil¬ 
ity p for a given settled system to colonize a particular 
neighboring stellar system (“sociological drive”), and the 
maximum travel distance D max that its spacecraft can 
travel, which thus determines adjacency (“technological 
prowess”). The Monte Carlo simulations of these per¬ 
colation processes allowed the calculation of the number 
of systems originating colonization programs that could 
reach the Solar System as a function (j>{p,D max ). In 
other words, <f(p, D max ) is a correlation function reflect¬ 
ing whether the Solar System and any other given star 
system are in the same cluster of systems colonized by a 
given interstellar civilization. 

This paper asks a different question, namely: for this 
calculated function <f>(p, D max ) and the evident lack of 
visitors to the Solar System, what are the upper limits 
on the likelihood of spacefaring life originating on a given 
habitable planet, as a function ofp and D max ? The num¬ 
ber of interstellar civilizations reaching the Solar System 
must be the number of stellar origin sites for colonization 
efforts could possibly reach us, multiplied by the proba¬ 
bility of a star-faring civilization arising at each star sys¬ 
tem. From the self-evident lack of visitors - meaning the 
number of colonizing efforts reaching the Solar System is 
less than one - the probability of a species willing and 
able to execute a colonization program must be simply 
the inverse of the probability of such civilizations evolv¬ 
ing in a given star system. Section |TT] develops this idea 
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further, by using an equation for the expected number 
of interstellar civilizations based on a break-down of the 
various relevant factors in a manner similar to the Drake 
equation, and solving for the probability of a spacefar- 
ing species arising on a given habitable world in terms of 
current constraints on such exoplanets. As with Ref. 4, 
all types of stellar systems are lumped into three broad 
categories - bright single stars, dim single stars, and mul¬ 
tiple star systems. Specific forms of 4> a (p, D max ) are ob¬ 
tained by fitting the results of the simulations in Ref. 4 
to a modified power law, where a runs over the three 
system types. Using these fitting functions, an example 
case is shown for a choice of habitable planet occurrence 
rates around various star system types, and the maxi¬ 
mum possible probability of a spacefaring species arising 
within the local Solar neighborhood are calculated as a 
function of p and D max . These results are summarized 
and discussed in Section ITTT1 


II. CONSTRAINTS ON LOCAL 
INTERSTELLAR CIVILIZATIONS 


As mentioned in Section |TJ this work begins with 
an equation based on the Drake equation, which 
parametrizes the number of extraterrestrial civilizations 
in the Milky Way galaxy from a combination of astro- 
physical, biological, technological and sociological inputs, 
with an eye towards estimating the total number of civi¬ 
lizations which are available for radio or other communi¬ 
cations. A similar idea is useful in estimating the possi¬ 
bility of direct contact with an alien civilization. Specif¬ 
ically, such visits are only possible if (1) a potentially 
space-faring civilization arises in another star system, 
and (2) this civilization has the willingness and technol¬ 
ogy to reach the Solar System, possibly after colonizing 
one or more star systems along the way. It is now pos¬ 
sible to make educated guesses on the likelihood of the 
first point. The result of this computation is a variable 
Nk £ [0,1], for a star system k, where denotes the 
probability of an interstellar civilization arising in star 
system k. Note that we use the term “system” here, since 
we include the possibility of space-faring species origi¬ 
nating within a multiple star system. The second is an 
appropriate setting for the percolation model developed 
earlier in Ref. 4. The results calculated in that paper 
implicitly give a correlation function between the Solar 
System and all other systems in question, e.g. the prob¬ 
ability that the Solar System and another given system 
are both within the same colonization cluster, based on a 
choice for colonization probability p and travel distance 
Dmax- For a stellar system k , this correlation probability 
is denoted <j>k(p, D max ). Then the expected number (N) 


of civilizations visiting the Solar System is given bid 

Nkfik {jPi Dmax) ( 1 ) 

k 

We can parametrize Nk much like the Drake equation to 
get 

w = E f P ,knh,kfe,kfs,k^k(p , D max ) (2) 

k 

where the parameters f p ,rih,fe, and fs, for each stellar 
system k, are defined in Table [T] 

Symbol Definition 

f p Fraction of stars with planets 
nh Number of those planets within 

the habitable zone of their stellar system 
fe Fraction of planets in habitable zone that 
develop life 

fs Fraction of life-bearing planets developing 
a star-faring civilization 

TABLE I. Definitions of parameters used in equation § for 
the expected number (N) of interstellar civilizations reaching 
the Solar System. The total probability Nk of a space-faring 
civilization arising in a given star system is the product Nk = 
fpUhfifs- 

Since Ref. 4 only considers those star systems within 
40 pc, we limit our discussion to that volume to find the 
expected number of civilizations reaching the Solar Sys¬ 
tem. Theoretically, at least, the astrophysical properties 
of each stellar system within that volume - i.e. the num¬ 
ber of planets in each system, and the number of those 
within the habitable zone for that system - could be di¬ 
rectly measured, or at least estimated based on the age, 
spectral type, and metallicity of each star present in the 
system TO]. However, even those projects such as the 
Kepler spacecraft have garnered a huge amount of infor¬ 
mation about exoplanets, the study of other worlds is 
still young, and only broad outlines can be seen at the 
moment. Thus, to match current knowledge, the sum 
© for the expected number of visiting civilizations is 
divided into three pieces, according to whether each sys¬ 
tem consists of a single star of spectral class K or earlier 
(“bright”), a single star of class M or other non-stellar 
object (“dim”), or a system of multiple stars (“multi¬ 
ple”). 

In addition, there are two parameters more in the as- 
trobiological realm than the astronomical, namely fa and 
fs- The fraction fa of habitable zone worlds where life 


1 Equations (2) and (3) of Brin [2] are similar in spirit to the 
equations presented here, but use different parameters to define 
the model space. 
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evolves is unknown at the momen© although this could 
be constrained by potential near-term space missions, 
such as those craft currently exploring Mars, or possi¬ 
ble missions for Europa, Titan and Enceladus. On the 
other hand, the fraction fg of inhabited worlds where 
an interstellar civilization arises is completely unknown, 
since as far as we know, the answer is zero. In light of 
our ignorance of the processes behind fe and fg, it is 
reasonable to keep the model as simple as possible; thus, 
for the following the product fefs is assumed to be the 
same for all star system types, and the remaining param¬ 
eters are taken to be identical for a given system type 
a £ {b,d,m}. This means the sum ([2]) for the expected 
number of visiting civilizations is given by 


(N) = (Ms) 


^2(fpn h ) a ^2(l)k(p,D 

r 


k£a. 


The results of Ref. 4 allows the functions 

^ ^ 4*k{.Pi D max ) 


to be calculated for all systems k within 40 pc of the Solar 
System with a particular system type a , so the previous 
equation Q for (N) will be used in the form 

(N) = (fefs) (fp n h) a (N a (p,D max )) (3) 

ae{6,d,m} 



FIG. 1. Contour plot of the maximum possible probabil¬ 
ity fefs of a star-faring civilization arising on an Earth-like 
planet in another star system, with spacecraft from that civ¬ 
ilization reaching the Solar System. This graph uses the 
choice of Eartli-like planets existing in other star systems 
(fpTih)a = (0.0232,0.309,0.0174), for bright and dim stars, 
and multiple star systems, respectively. 


As a way of seeing the behavior of the number of sys¬ 
tems that can reach the Solar System, we fit each of the 
three functions (N a ) to a modified power law in p, given 
by 


This equation for the expected number of interstellar 
civilizations to reach the Solar System is not helpful as it 
stands, since it depends on unknown probabilities fe and 
fg. However, the only observable fact that can be used 
in this discussion namely, the complete lack of evidence 
that the Solar System has been visited by extraterrestrial 
spacecraft - does serve to constrain the product fefs- In 
other words, the fact that (N) < 1 allows the calculation 
of upper limits on the product fefs, since this limit on 
the maximum value of the equation © for (N) implies 


(Ms) < 


'y ' (fpHh)a(N a (p, Drnax)) 


L aG{fc,d,m} 


-1 


( 4 ) 


(N a (p, D max )) = A aP c ^ D ™*+( 5 ) 

This gives us the coefficients given in Table [TTJ where the 
values of (N a ) are fit for p € [0,1] and D max € [3.0, 7.0] 
pc. The large and negative coefficient of the linear D max 
terms in the exponent show that, as the maximum travel 
distance increases, a smaller value of the colonization 
probability is necessary in order to achieve the same over¬ 
all value of (N a ), for a given system type. However, the 
positive quadratic coefficient shows that this tendency 
levels off a bit as D max increases. Both of these effects 
are more pronounced for the bright stars and multiple 
star systems than they are for the dim stars. 


All quantities on the right-hand side of this inequality 
are either probabilities now grounded in exoplanet sur¬ 
veys (f p ) and climate models (tih) for these worlds, or 
else calculated from the percolation model ((N a )) of Ref. 
4, based on assumed values for the sociological drive and 
technological abilities of an interstellar civilization, i.e. 
the colonization probability p and maximum travel dis¬ 
tance D max , respectively. 


2 As noted by Spiegel and Turner [9], current knowledge about 
the evolution of life (specifically on Earth) does little to limit the 
value of ft. 


System type 

A 

ci (pc -2 ) 

c 2 (pc -1 ) 

C3 

bright 

2074 

0.4905 

-5.743 

18.52 

dim 

1598 

0.2351 

-2.714 

9.668 

multiple 

1575 

0.4896 

-5.753 

18.60 


TABLE II. Fitting coefficients for the functions 
(N a (p, Dmax)) giving the number of stellar systems of 
a given type a £ { b , d, m} able to reach the Solar System for 
a given choice of colonization probability p and maximum 
travel distance D max (in pc). The functional form of 
(N a (p, Dmax)) is given in equation 0. 

As a specific instance of the use of this model, values 
for the parameters are chosen based on a geometric mean 
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of results given in the literature, cited by Winn and Fab- 
rycky [TO] , Specifically, the occurrence rates from Table 
2 of that work are used, where the three listings for M 
stars are used for dim stars, the four entries for FGK 
are combined with the GK entry for bright stars, and 
the occurrence rate for habitable planets in multiple star 
systems is chosen as 75% of that for bright stars. This 
gives the parameter choices of 

{fp^h)(y = l(fp'Tlh) bright > ^fp^h) dimi {fp'hlh) multiple] 

= (0.0232,0.309,0.0174) 

The maximum possible values of fifs are then calcu¬ 
lated, and a contour plot showing the results is given 
in Figure [T] We now discuss the results shown in this 
figure; these comments relate specifically to the choice 
of ( f p nh)a used here, but the broad outline remains the 
same regardless of what values are chosen. 

• Interstellar civilizations capable of making longer 
journeys between star systems (higher D max ), and 
more likely to begin these journeys (larger p ), have 
a higher chance of reaching the Solar System. Thus, 
the empirical fact that visiting spacecraft of this 
type are not observed places stronger limits on 
fifs , and so the upper limit of this product de¬ 
creases in size as we move towards the upper-right 
hand corner of Figure |Tj 

• Specifically, the maximum range of travel distances 
for interstellar colonizers studied in this research 
was chosen to be 7.0 pc. If such a civilization al¬ 
ways colonized its neighbors (p = 1), then the up¬ 
per limit on the probability of such a species arising 
within a biosphere less than 40 pc away from the 
Solar System is fifs < 0.00171. In other words, at 
most only one of out every 585 habitable planets 
within the local Solar neighborhood could be the 
cradle of an interstellar civilization in order not to 
have evidence of their spacecraft in our Solar Sys¬ 
tem. 

• There is a relatively large regime of colonization 
probabilities p and maximum travel distances D max 
where interstellar civilizations would never reach 
the Solar System, even if technological civilizations 
capable of star travel arose on every single Earth- 
like planet in the Solar neighborhood. The reach 
of these interstellar civilizations could be quite siz¬ 
able, as seen in Ref. 4. However, the practical ef¬ 
fect here is that there is a region in p — D max space 
where ( N ) < 1 regardless of the value of fifs- 

Recall that the product fifs deals with the chances 
of an interstellar civilization arising on a given habitable 
planet, where we have estimated values for the expected 
number of such planets per stellar system (as divided be¬ 
tween bright and dim stars, and multiple systems). Since 
the latter is on the order of 1-10%, then the most strin¬ 
gent limit on the probability of a space-faring species 


evolving within any star system is around 10 -5 . Rela¬ 
tive changes in the sizes of (f p rih)a for the three different 
system types a may change the shapes of the contours 
given in Figure [l] but the upper limit values for a given p 
and D max will be altered only for absolute changes in the 
order of magnitude for the number of habitable planets 
per system, regardless of system type. 


III. CONCLUSIONS 

In this paper, the model of interstellar colonization as 
a percolation process, originally developed by Landis and 
applied to the physical positions of all star systems within 
40 pc, is used to develop upper limits on the probability of 
an interstellar civilization arising on an Earth-like planet, 
and spacecraft of that civilization reaching the Solar Sys¬ 
tem. The model was parametrized in a manner similar to 
the Drake equation, with coefficients giving the fraction 
f p of stars with planets, and the number of those plan¬ 
ets rih within the habitable zone of each stellar system, 
the fraction of habitable planets fi that actually develop 
life, and finally the fraction of those biospheres fs which 
originate a space-faring civilization. At this point in our 
knowledge, there are now empirical constraints on the 
first two parameters, but the latter are more or less free. 
The philosophy of this paper is to use a simple model in 
order to place some upper limits on those values dealing 
with the origin of life and interstellar civilizations, since 
there is currently no observational data on these issues. 
Our ignorance of many of the processes involved in the 
evolution of a technological species, and the likelihood of 
that species choosing to develop interstellar settlements 
leads to the choice of the simple model given here, with 
as few parameters as possible to give a reasonable look 
at the question at hand. 

However, in some sense we have just moved our igno¬ 
rance - instead of leaving ft and fs open, this work puts 
upper limits on those values at the expense of two new 
coefficients, the probability p of a space-faring species 
to colonize one of its neighboring star systems, and the 
maximum distance D max that can be considered ” neigh¬ 
boring” . Yet this shift in focus may be helpful, in that it 
puts a different light on questions relating to the Fermi 
paradox. Indeed, as shown elsewhere |3j, a maximum 
travel distance D max of about 3 pc is sufficient to allow 
effective interstellar travel, i.e. actual travel paths are 
less than twice the shortest distance between end points, 
with the benefit of passing through intermediate star sys¬ 
tems, allowing for replenishment and refit. As for the 
colonization probability p, as yet humanity itself has not 
taken this step, but it is for each of us to evaluate the 
likelihood of this happening when technology has reached 
the point where this question is a viable one. 

The percolation model presented here, however, only 
includes nodes representing all star systems with 40 pc of 
the Solar System, and thus does not account for the pos¬ 
sibility of interstellar civilizations arising outside of that 
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volume. Indeed, it is perfectly conceivable that a galaxy¬ 
wide colonization effort may have started from any point 
within the Milky Way, as long as the colonization prob¬ 
ability p is large enough, leading to much smaller upper 
bounds on fefs- To truly study this issue, the perco¬ 
lation model would also have to be extended onto the 
galactic scale. The stellar number density - the expected 
number of star systems found inside a cubic parsec of vol¬ 
ume, for example - is not constant over the Milky Way, 
but varies according to both distance away from the cen¬ 
ter of the galaxy, as well as scale height away from the 
galactic plane IT]. Thus, the study of percolation clus¬ 
ters using a physically realistic distribution of simulated 
stellar systems would not correspond to previous mod¬ 
els from statistical mechanics. Said in a different way, 
the critical probability p c , i.e. the probability to have a 
galaxy-spanning cluster, needed as a function of D max 
is currently unknown. Thus, the size of the largest clus¬ 
ter, corresponding to the interstellar civilization with the 
greatest reach, as a function of the colonization proba¬ 
bility p and the maximum travel distance D max between 


two adjacent systems, remains a problem to be studied 
in future simulations. However, having a relatively low 
value of p c would imply that the upper limits on fofg 
would be very small, or else a galactic civilization would 
have arisen and easily reached out Solar System, contrary 
to observation. 

Another important caveat for the work here is that 
none of the relative motion of star systems in the Solar 
neighborhood is taken into account. This motion leads 
to changes in positioning over periods of only 10 1 2 3 — 10 4 5 6 
yr The Monte Carlo simulations used to model the 
percolation process in Ref. 4 are “timeless”, in the sense 
that all systems are static, and the speed of colonizing 
spacecraft is not considered. It is likely this time scale 
is roughly the same as the process for a newly settled 
star system from an interstellar civilization to built its 
own colonizing spacecraft, and then this craft reaching a 
new system, so the dynamical changing of available star 
systems for a star-faring species would definitely affect 
the results presented here; this also is an issue for further 
work. 


[1] Bailer-Jones, C. A. L. 2014, accepted in Astronomy & 
Astrophysics, arXiv: 1412.3648 [astro-ph]. 

[2] Brin, G. D 1983, Quarterly Journal of the Royal Astro¬ 
nomical Society, 24, pp. 283-309. 

[3] Cartin, D. 2010, Journal of the British Interplanetary So¬ 
ciety, 63, pp. 218-221, arXiv: 1104.4012 [physics-pop.ph]. 

[4] Cartin, D. 2014, Journal of the British Interplanetary 
Society, 67, pp. 119-126, arXiv:1404.0204 [physics.pop- 
ph[. 

[5] Hair, T. W., Hedman, A. D. 2012, International Journal 
of Astrobiology, 12, pp. 45-52. 

[6] Hart, M. 1975, Quarterly Journal of the Royal Astro¬ 


nomical Society, 16, pp. 128-135. 

[7] Juric, M., Ivezic, Z., Brooks, A. et al. 2008, Astrophysical 
Journal, 673, pp. 864-914, arXiv: 0510520v2 [astro-ph], 

[8] Landis, G. A. 1998, Journal of the British Interplanetary 
Society, 51, pp. 163-166. 

[9] Spiegel, D. S., Turner, E. L. 2011, Proceedings of the 
National Academy of Sciences, 109, pp. 395-400, arXiv: 
1107.3835v4 [astro-ph.EP]. 

[10] Winn, J. N., Fabrycky, D. C. 2015, Annual Reviews of 
Astronomy and Astrophysics, 53, arXiv: 1410.4199 [astro- 
ph.EP], 



